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With a treatment of the 4/ states of EuB6 based on LDA+U method, the mixing of Eu / states 
with B p states around the X point of the Brillouin zone is shown to have unexpected consequences 
for the effective exchange interactions. We analyze in detail the orbital character of electronic states 
close to the Fermi level and discuss the effective exchange between the itinerant electrons and the 
local 4/ moments. The analysis suggests that the ordered phase may provide the first example 
of a half metallic semimetal, and that the physics of EuB6 should be described in terms of a two 
band Kondo lattice model with parallel (ferromagnetic) coupling of the conduction electrons and 
antiparallel (antiferromagnetic) coupling of the valence electrons to the local 4/ moments. 
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I. INTRODUCTION 

Despite numerous experimental and theoretical stud- 
ies, understanding of the electronic properties and mag- 
netic coupling of EuBg still provides challenges. While 
having a simple crystal structure, consisting of a simple 
cubic lattice of Eu atoms with a Bg octahedron located 
in the center of each cubic cell, nevertheless the trans- 
port and magnetic properties of this system are complex. 
EuBg orders ferromagnetically at 15.1 K, which is accom- 
panied by a huge decrease of resistivity and a significant 
blue shift of the reflectivity plasma edge. 1 At 12.7 K an- 
other phase transition takes place, which is observed as 
a broad peak in the specific heat or an anomaly in the 
resistivity^ The origin of this transition is still unclear, 
with possible explanations including spin reorientation 3 
or a long-wavelength modulation of the spin density. 4 
Besides these properties, EuBg exhibits a rather sluggish 
increase of magnetization with decreasing temperature 5 
and unusual pressure dependence of the Curie tempera- 
ture with strong increase up to 70 kbar and a flat depen- 
dence at higher pressures £ 

The LDA electronic structure of EuBg was previously 
investigated by Hasegawa and Yanase 7 - and Massida et 
al. 8 Small overlap of conduction and valence bands was 
found resulting in Fermi surface pockets centered at the X 
point of the Brillouin zone. Isostructural CaBg and SrBg 
have similar bandstructures, which have seemed to be 
consistent with their observed transport properties^ (al- 
though they are not understood in any detail). Recently 
CaBg has been studied with several variants of the GW 
method: the conventional pseudopotential GW (iS and 
all-electron GWjii and a self-consistent GW method. 12 
While there are differing results among these, the ma- 
jority of them predict the opening of a gap of the or- 
der of 1 eV, making the LDA conclusions about the 
groundstate of divalent hexaborides questionable. Re- 
cent angle-resolved photoemission measurements^*^ re- 
ported a bandgap in CaBg, and only recently has it 
been demonstrated that synthesis from ultrapure boronA 5 " 



leads to transport properties that are characteristic of a 
semiconductor rather than a semimetal. For the param- 
agnetic phase of EuBg, a bandgap was also observed^ 
below an occupied electron pocket at the X point that 
was interpreted as carriers resulting from (a high density 
of B) vacancies. Recently the transport properties have 
been interpreted in terms of a similar models 

As usual in LDA-based band structure calculations of 
rare earth systems, the 4/ states must be treated in a spe- 
cial way in order to insure their correct filling. In order to 
enforce the correct filling the 4/ states were treated sep- 
arately from the rest of the system in Ref. H, neglecting 
any hybridization involving the / states and thus giving 
no insight into magnetic coupling mechanisms in EuBg. 
Here we use the LDA+U method, which obtains the cor- 
rect 4/ state filling while keeping these 4/ states in the 
same Hilbert space as the rest of the system, and so al- 
lows for mixing of 4/ states with the valence states. The 
aim of the present work is to investigate the effective ex- 
change interaction between the localized / moments and 
the band electrons. For this purpose we perform a de- 
tailed analysis of the orbital character and dispersion of 
the states close to the Fermi level. Based on this analysis 
we suggest an unusual two band Kondo lattice model to 
be the relevant picture for understanding the magnetic 
behavior of EuBg. 

II. COMPUTATIONAL DETAILS 

The calculations were performed using the full- 
potential linearized augmented-plane-waves (FLAPW) 
method as implemented in the Wien2k code>i£ The 
LDA+U method was used with the double-counting 
scheme of Anisimov and collaborators^ The standard 
parametrization of the on-site Coulomb interaction in- 
volves two parameters U and J. However, since we are 
dealing with the Eu 2+ / 7 ion, which has completely filled 
spin-up / shell and completely empty spin-down / shell, 
the role of J reduces to merely renormalizing the U value. 
Therefore we can set J=0 and quote just the U value. 
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All the presented calculations were performed with LDA 
exchange-correlation potential in the parameterization of 
Perdew and Wang. -4 The calculations were performed 
without spin-orbit coupling in the scalar relativistic ap- 
proximation as implemented of Wien2k code. 

The following computational parameters were used. 
The atomic radii were 2.7 and 1.5 Bohr for Eu and B re- 
spectively. The APW+lo basis se*2S, with additional lo- 
cal orbitals for Eu 5s and 5p states, was characterized by 
plane-wave cut-off R m tK max — 7. We used 56 k-points 
in 1/48 irreducible wedge of the Brillouin zone. The nu- 
merical convergence of the total energy was better than 
0.1 mRy. The internal parameter was relaxed so that the 
corresponding force was smaller than 1 mRy/a.u. 
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FIG. 1: The total energy vs. volume curve calculated 
III. RESULTS AND ANALYSIS with U=7eV. In the inset the relative change of the nearest- 

neighbor (dotted) and intra-octahedron (solid) B-B distance 
A. Bulk properties as a function of relative change of the lattice constant. 



We have performed total energy vs. volume calcula- 
tions in order to (i) determine the theoretical equilib- 
rium volume and bulk modulus and (ii) investigate the 
influence of pressure on the band structure. For each 
value of lattice constant we have optimized the nearest- 
neighbor B-B distance (the only free internal parameter) 
using atomic forces. The calculated energy vs. volume 
curve is shown in Fig. ^ The inset shows the scaling of 
B-B distances with the lattice constant indicating that 
the bonds within B-octahedra are more rigid than the 
nearest-neighbor (inter-octahedra) B-B bond. The calcu- 
lated equilibrium volume of 70.1 A 3 is about 96 % of the 
experimental valued The calculated bulk modulus of 161 
GPa is agrees very well with the experimental valued 
The present results were obtained with U=7 eV, how- 
ever similar calculations performed for U=6, 8, and 9 eV 
(without optimization of the internal parameter) showed 
that the bulk properties are insensitive to the value of U. 



B. Bandstructure and exchange coupling 




FIG. 2: Spin-polarized density band structure obtained with 
U=7eV at the experimental lattice constant. The spin-up 
bands are marked red (brighter). 



In Fig. |21 we show the spin-polarized bandstructure 
obtained at the experimental lattice constant with U=7 
eV. Taken literally, the bandstructure indicates a metallic 
ground-state with band overlap (negative gap) around X 
point in both spin channels. More importantly the over- 
lap of majority (spin up) bands is strongly enhanced in 
comparison to minority (spin down) bands. The origin 
of the different band overlap is the opposite sign of the 
up/down exchange splitting induced in the conduction 
and valence band. While lowering of the energy of the 
spin-up conduction band with respect to the spin-down 
band points to parallel coupling of conduction electrons 
to the local / moments, the energy of the spin-up valence 
band is higher than that of the spin-down band indicat- 
ing antiparallel coupling to the local / moments. In the 



following we investigate in detail origin of this particular 
effective exchange coupling. There are 6 symmetry re- 
lated X points in the Brillouin zone of simple cubic struc- 
ture. The following analysis is performed particularly for 
X=(0, 0,1/2) and all references to spatial orientation the 
orbitals are with respect to this choice. Also, we do not 
distinguish between directions in real and k-space, which 
is not confusing for orthogonal unit cell. 

The modulus of the wave-function corresponding to 
valence band maximum at the X point is shown in Fig. 
[3] In Fig. 0] we show the schematic plot of the same 
wave function in terms of B-p and Eu-/ orbitals. Let us 
make several observations: (i) p-orbitals of B atoms at 
(1/2, 1/2, z) do not contribute to this wavefunction while 
the other B atoms (which lie in the z—1/2 plane) con- 
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FIG. 3: Contour plot of the absolute value of the minority- 
spin valence band wavefunction at the X point. The right 
panel shows a cut through the boron plane (001) perpendic- 
ular to T-X direction. The left panel shows a cut by (110) 
plane going through the center of the unit cell. Mixing with 
/ state in corners of the plot is visible even for the minority 
spin. This feature is much stronger for the majority spin. 




FIG. 4: Schematic plot of the valence band state at the X 
point showing the phases of participating orbitals. Two unit 
cells in X-M direction are shown. 



tribute their p x or p y orbitals, (ii) hybridization between 
p orbitals in the neighboring unit cells (in X-M direction) 
is forbidden by symmetry at the X point, (iii) dispersion 
in the T-X direction in enhanced by hybridization with 
Eu / orbital of xyz symmetry (note that this hybridiza- 
tion is allowed due to a phase shift of tt between adjacent 
Eu layers and thus is forbidden at T point). Together 
points (ii) and (iii) explain the convex dispersion of the 
valence band with the top at the X point. The hybridiza- 
tion with Eu / also explains the antiparallel coupling of 
the valence electrons to the local / moments. The spin 
up valence orbital hybridizes strongly with the occupied 
/ state localized close to the Fermi level and thus its en- 
ergy is increased due to band repulsion. The hybridiza- 
tion shift (level repulsion) in the spin down channel is 
much weaker and of opposite sign, since the unoccupied / 
are localized high above the Fermi level due to the on-site 
Coulomb repulsion. As a result we obtain an effective an- 
tiferromagnetic exchange interaction of kinematic origin 
in the way described by the periodic Anderson models 
Similar analysis is performed now for conduction band. 
By analyzing the orbital contributions to the conduction 
band we find that it contains a mixture of B-p and Eu-rf 
states, where the d content decreases when going away 
from the X point and vanishes at the T point. We show 




FIG. 5: Contour plot of the absolute value of the minority- 
spin conduction band wavefunction at the X point. The left 
panel shows a cut through the Eu plane (001) perpendicular 
to T-X direction. The right panel shows a cut through the 
center of the unit cell by the (100) plane. 




FIG. 6: Schematic plot of the conduction band state at the 
X point showing the phases of participating orbitals. The left 
panel shows the view along the F-X direction. The right panel 
shows the view along the X-M direction. 



the calculated modulus of the wave-function correspond- 
ing to the conduction band at the X point in Fig. Eland 
the corresponding schematic plot in Fig. We make fol- 
lowing observations: (i) the only orbitals that contribute 
significantly are p z orbitals localized on B atoms in the 
z=l/2 plane and Eu d x 2_ y 2 orbitals, (ii) at the X point 
the p orbitals in the neighboring unit cells (in X-M direc- 
tion) form a bonding orbital, which can hybridize with 
Eu d x 2_ y 2 orbital, going away in the X-M direction in- 
troduces a phase shift between p orbitals in neighboring 
cells a reduces hybridization with the d states, eventually 
at the M point the p orbitals form an antibonding com- 
bination and hybridization with d's is suppressed, (iii) 
looking along X-r direction we find that the phase shift 
of 7T at X point allows the p states to hybridize with Eu 
d x 2_ y 2 states while this mixing is forbidden at T point. 
In conclusion the points (ii) and (iii) explain the concave 
dispersion of the conduction band with bottom at the 
X point. The conduction band originates from the B p 
band which mixes strongly with Eu d x 2_ y 2 band close to 
the X point. This picture is corroborated by the band- 
structure of the empty boron lattice (without Eu atoms) 
shown in Fig. \7\ which contains a similar conduction 
band but with greatly reduced dispersion. The exchange 
interaction of the electrons in the conduction band with 
local / moments is of ferromagnetic / — d intra-atomic 
origin. Both the reduced value of the exchange splitting 



4 




FIG. 7: The band structure of empty boron cage. The correct 
filling (2 additional electrons per octahedron) was achieved by 
virtual crystal approach. 

as compared to higher lying d-bands and decrease of the 
exchange splitting when going away from X point are 
easily explained by variable content of the Eu d in the 
conduction band states as described in (ii) and (iii). 

C. Role of Pressure and U 

In the previous section we have shown that the mech- 
anism of exchange with the local / moments is rather 
different for conduction and valence electrons. The con- 
duction electrons are polarized through intra-atomic ex- 
change interaction between the Eu d and / states. The 
strength of the interaction is determined mostly by the 
content of the d orbital in a particular wavefunction and 
is independent of the energy of / states. The polariza- 
tion of valence electrons arises from different hybridiza- 
tion splitting in spin-up and spin-down channel. Since 
the hybridization is negligible in the minority channel, 
the strength of the effective exchange is determined by 
hybridization shift of the majority valence states. This 
shift is inversely proportional to the energy difference be- 
tween the / and valence band. The position on the energy 
scale of the / states is therefore crucial for determining 
the strength of the effective exchange. 

The position of the lower and upper 4/ bands is de- 
termined by two factors within the LDA+U approach: 
(i) the center of gravity of the LDA / bands, (ii) the 
screened on-site Coulomb interaction U. The center of 
gravity of the / bands depends on charge transfer, but 
its LDA position suffers from well-known self-interaction 
error. The additional terms in the LDA+U hamiltonian 
enforce the splitting into lower and upper Hubbard band 
and in an approximate way correct for the self-interaction 
error in the lower band. To calculate the precise value of 
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FIG. 8: The band overlap of the spin-up bands as a function 
of the parameter U evaluated at the experimental volume is 
shown in the left panel. The band overlap in the spin-up 
channel (upper curve) and spin-down channel (lower curve) 
as a function of unit cell volume evaluated at U=7 eV. 

U, besides the fact that it is defined in a somewhat loose 
sense, is beyond the scope of this work. The typical U 
quoted for 4/ electrons after accounting for screening is 
between 7 and 9 eV. In Fig. |H|we show the band overlap 
in the majority spin channel as a function of U from 6 
to 9 eV. We point out that U around 7 eV yields the / 
bands approximately 1 eV below the Fermi level, which 
is the position deduced from optical experiments^ and 
therefore we assume the results for this value to be most 
realistic. 

The effect of applied pressure amounts to overall 
broadening of the bands and thus enhanced overlap of 
the conduction and valence bands at the X point. For a 
unit cell volume increase of about 10% a gap in the mi- 
nority spin channel opens, while an appreciable overlap 
remains in the majority spin channel (see Fig. [SJl. 

D. Model Hamiltonian 

In the previous sections we have described the origin 
of dispersion and effective exchange coupling with local 
/ moment in both valence and conduction bands. Based 
on the band structure we suggest the following model 
hamiltonian to provide a reasonable description of the 
low energy physics: 

H = J2 K v l<? v *° + 44.^] - 

k,<r 

1 l J iMk' S i ■ vl a (T a pv k ^+ (1) 

i,kk' 

^kk' S i ' cl a eT a pc k 'i3], 

where operators Vk a and Cko- correspond to valence and 
conduction bands respectively, Si is a total spin operator 
of the local / moments and L is the number of unit cells 
in the normalization volume. 

Obtaining parameters of the model hamiltonian from 
the bandstructure is not a straightforward task, in par- 
ticular the choice of the reference non-interacting state is 



5 



tricky. It was shown by Schiller and Nolting 24 and fur- 
ther discussed by Miiller and NoltingS^ that in case of 
k-independent exchange parameter J, the majority-spin 
band is rigidly shifted with respect to its non-interacting 
counterpart, while the minority-spin band is modified be- 
yond the rigid shift due to interaction with the local mo- 
ments. The situation in EuBg is more complicated since, 
as discussed below in detail, the exchange parameter is 
strongly k-dependent for purely chemical reasons (hy- 
bridization). It is therefore not possible to distinguish 
unambiguously the deviation of the exchange splitting 
A(k) from a rigid band shift arising from the correlation 
effects in the minority-spin band from that originating in 
the 'chemical' k-dependence of the exchange parameter 
J. Since the k-dependence of the exchange splitting A(k) 
along the X~r and X-M directions is very pronounced and 
can be qualitatively understood in terms of k-dependent 
exchange parameter J we neglect the correlation effects 
in the minority spin channel. 

The non-interacting dispersion relations e k and e k ob- 
tained from the bandstructure (Fig. are considered 
separately. The interaction of the conduction electrons 
with the local moments is dominated by the / — d intra 
atomic exchange depending only on the d content in a 
particular conduction state. The 'non-interacting' con- 
duction band is then approximated by an average of the 
spin-up and spin-down bands. Since the hybridization in 
the majority band is the source of the exchange splitting 
in the valence band, the 'non-interacting' valence band 
is well approximated by the spin-down band. The Fermi 
level is located close to the top, resp. bottom, of the 
valence, resp. conduction, band and so the low energy 
bandstructure can be parametrized by anisotropic effec- 
tive masses and the band overlap. The effective mass 
tensor for an ellipsoid of revolution is characterized by 
only two independent parameters fi± and correspond- 
ing to X-M and X~r dispersion respectively. The effec- 
tive masses obtained by parabolic fit from the bandstruc- 
ture (Fig. EJ) are fi c ± = 0.23, /ijj = 0.47, (i" ± = 0.25, 
and = 2.2. The band overlap of the 'non-interacting' 
bands at X point is 0.34 eV. We point out that using the 
spin-up conduction band as a non-interacting reference 
(rigidly shifted) does not lead to significant modification 
of the effective masses. 

Since the mechanisms of effective exchange with local 
moments are different we have to discuss determination 
of the corresponding coupling constants J^kk' separately. 
The kk' dependent coupling constants can not be deter- 
mined directly from the bandstructure and additional as- 
sumptions must be made. We start with the conduction 
band. 

We mentioned earlier that the conduction band cou- 
pling to the local moment is mostly due to intra-atomic 
f — d exchange described by 

H df = - J ^ Sj ■ d\ a cr aS3 d iSil (2) 
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FIG. 9: The k-dependent exchange splitting: splitting of the 
conduction band is shown in the upper panel, the splitting 
of the valence band is shown in the bottom panel. The ex- 
change splitting calculated with the tight-binding expression 
for V/p=24 meV is marked with the dotted line. 

which leads to 

•W = Jc i(k - k ' )R '. (3) 

However, the conduction band is not a pure d-band, but 
contains a mixture of d and p orbitals Ck = a(k)c?k + 
(3(k)pk. Treating Hdf as a first order perturbation leads 
to an effective coupling of the conduction of the form 

JJikk/ = Je I ( k - k ') R '«(k)a(k'). (4) 

Here a can be chosen at convenience as real and positive. 
Finally we make connection to the band structure by 
observation that the k-dependent exchange splitting of 
the conduction band is given by 

A c (k) = e^-4 T = 2J5« 2 (k). (5) 

Finally putting Q and JSJ) together we obtain 

25J^ kk , = e^-k')** ^A c (k)A c (k'). (6) 

The k-dependent exchange splitting of the conduction 
band is shown in Fig. Vanishing of the exchange 
splitting when the T and M points are approached can 
be understood in terms of simple nearest neighbor (NN) 
hopping picture. The k-dependent hybridization between 
p and d bands vanishes at T and M, and thus the valence 
states have a pure p character at these points. The p— d 
NN hopping in the conduction band has the same sym- 
metry as the p — f hopping in the valence band leading 
to the same k-dependence of the hybridization 10 and 
(frU|l discussed below. 
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For the valence band the effective exchange coupling 
arises by Schrieffer- Wolff transformation^ from Ander- 
son lattice hamiltonian with the f — p hopping term 



H pf = ^2 ^ik//wk + h.c. 

i,k,(7 



(7) 



The corresponding exchange parameter for states close 
to the Fermi level (i.e. far enough from the / band) is 
given by 



JV 

J i.kk' 



~Vi ti V ik f[—^— + 

2 L £k - £/ £k' - £/ 



1 



f 



(8) 



£k - £/ - U e k < - e f - U ] 



where e/ is the energy of the occupied / states and e/ + U 
is the energy of the unoccupied / states. As in the case of 
the conduction band we have to make additional assump- 
tion in order to extract information about J" kk / from the 
bandstructure. The natural assumption in terms of the 
local orbitals is to expect non-zero hopping Vf p only be- 
tween the nearest neighbor p and / (xyz) orbitals (see 
Fig. 0J. Note that the sign of Vf p , given by symme- 
try of the orbitals, is different for different pairs of or- 
bitals (each p orbital has 8 NN / orbitals). In particular, 
pairs connected by a vector (m,n, 1/2) have the same 
sign, which is opposite to the sign for pairs connected by 
(m, 7i, — 1/2) (m,n = ±1/2). Now we can write down the 
matrix element Vki as 



where 





= 2^e" 


ik-R; 




F(k) = 


sin(^ 4 


aky 

' ~Y 


+ a^ ) + 




sm(— 4 




_a_^ ) + 




' ~2~ 






sm(— - 


aky 

~Y 


+ ^) + 




. / cik x 

sm (— - 


' ~Y 





(9) 



(10) 



Providing we know the value of the hopping parameter 
Vfp the coupling J? kk , can be calculated by feeding the 
above expression and the known dispersion relation into 
equation JSJ. In order to obtain the parameter Vt p we 
express the exchange splitting of the valence band as 



A"fk) 



6 kT _t k| 



8SVf p F(k) 2 [ 



1 



£k - e f 



e k -tf-U 
(11) 

and use the fact that along the T-X direction F(k) = 
4sin(afc 2 /2), while along the X-M direction F(k) = 
Acos(ak x /2). The exchange splitting of 0.26 eV and 
e k — ej w 1 eV at the X point yield Vf p =24 meV. In Fig. 
[3] we show the k-dependent exchange splitting together 
with the expected k-dependence. For the X-M direction 



we show only the region close to the Fermi level. The 
reason is that close to the crossing of the non-interacting 
band with the / level the exchange splitting A(k) is not 
well defined. 



E. Mean field results 

Here we want to illustrate the potential usefulness of 
the suggested hamiltonian by investigating some of its fi- 
nite temperature properties in mean field approximation. 
Due to the three dimensionality and the large moment, 
mean field should be realistic in many respects. A simi- 
lar approach was taken by Korenblit 4 without providing 
detailed results. We focus on the role of the exchange en- 
hanced band overlap and its consequences. In the mean 
field approximation the hamiltonian Q reduces to 

Hmf - Yl [ £ k>L«k CT 4 eLcLc kCT ] - hS i> ( 12 ) 



where the spin-dependent dispersion relations are given 



by 



ka 



e£ - crJ a (S z ) and h = £ a J a (n^ - nf) is an 



effective self-consistent magnetic field, with a = v,c and 
a = ±1. The quantities n^n® are the occupations of 
the bands. Going to the hole picture of the valence band 
and using the charge neutrality condition to determine 
the Fermi level we obtain the following set of equations 

=0 (13) 
h = Y,<r{J°<(S*)- rnWj) (14) 

S Z =B V ^ (15) 

where B 7 / 2 {x) is the Brillouin function for S = 7/2. 

We have solved the mean field equations using the ef- 
fective masses listed above. To make the analysis as sim- 
ple as possible we approximate the coupling parameters 
J with a single k-independent constant (with opposite 
sign for the valence and conduction bands). To assess 
sensitivity of the model to the choice of J, we have used 
two values: (i) \J\ — 0.04 eV corresponding to the X 
point exchange splitting of 0.28 eV (see Fig. 2) and (ii) 
| J | = 0.10 eV used by Korenblit. 4 The calculated ordered 
moment vs. temperature S Z (T) curves are shown in Fig. 
8. There are two distinct regimes. 

(i) For band overlap A < (i.e. bandgap of |A|) there 
is a minimum ordered moment S ml „ necessary to estab- 
lish the overlap of spin-up bands. Below this value the 
magnetization cannot be self-stabilized since the corre- 
sponding effective field h is zero. Solution of equation 
(15) can be visualized as the intersection, on the inter- 
val (0,7/2], of the right-hand side s = B 7 / 2 as a func- 
tion of S z , given by (14), and the straight line s = S z 
of the left-hand side. The right-hand side is zero for 
S < S m i n and, from definition of -67/2(2;), less or equal 
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to 7/2. This means that curve representing the right- 
hand side must cross the line representing the left-hand 
side an even number of times since it is continuous and 
its initial and final points are in the same half-plane de- 
termined by the left-hand side. Therefore there must be 
an even number of non-zero solutions of (15) for A < 
and a given temperature. In our case it means either zero 
or two solutions. 

In order to identify, in the two solution case, the solu- 
tion with lower free energy we approximate the total free 
energy by the sum of the free energy of non-interacting 
valence and conduction electrons evaluated with the self- 
consistent value of the effective field h and the free energy 
of non-interacting local moments in the effective field h 
at a given temperature. In all cases we find that the 
solution with larger ordered moment corresponds to the 
lower free energy. 

The hamiltonian (1) preserves separately the num- 
ber of electrons in valence and conduction bands. If 
those bands are completely filled or empty in the non- 
interacting ground state, which is degenerate with re- 
spect to orientation of the local moments, there is no 
effective coupling between the local moments since the 
degenerate grounstate manifold is disconnected from the 
excited states. On the other hand when a minimum net 
magnetization already exists a band overlap is induced 
and a magnetic groundstate can be found as described 
by the mean field equations. There is no continuous con- 
nection between the non-magnetic and magnetic states 
and therefore we conclude that the transition is of the 
first order. 

(ii) For A > one solution exists in most of the stud- 
ied cases, with the ordering temperature increasing as 
A increases. The magnetization versus temperature be- 
havior deviates from the standard Weiss curve, which is 
obtained for linear dependence of the effective field h on 
S z . In a narrow range of A close to zero three solu- 
tions may exist, providing yet another phase transition 
below the magnetic-nonmagnetic one. We do not make 
any conclusions about the order of the phase transition 
in this parameter range. In the insets of Fig. 8 we show 
the ordering temperatures as a function of band over- 
laps. The difference between the ordering temperatures 
obtained with different coupling constants J indicates a 
strongly non-linear Tc vs. J dependence. (Note that if 
the band shifts due to the ordering of the local moments 
are very small, Tc oc J 2 Ref. 0). 



IV. DISCUSSION 

One of the controversial questions concerning divalent 
alkaline earth hcxaborides is whether the groundstate is 
insulating or metallic at stoichiometry. The experimen- 
tal as well as theoretical evidence is controversial. Re- 
cent LDA calculations by Massidda et a& neglecting the 
magnetic order yield a metallic band structure for EuBg 
similar to that of CaBg or SrBg with a small band over- 







in \ V 

y , iii. 




0.0 0.2 
\ A (eV) 

11 



2 4 
Temperature (K) 




10 20 30 

Temperature (K) 

FIG. 10: Ordered spin moment vs. temperature evaluated 
from the mean-field model for various values of the band over- 
lap. The upper and lower panel correspond to J=0.04 and 0.1 
eV respectively. The magnetization curves, ordered from the 
left to the right, correspond to band overlaps from -0.15 eV 
(bandgap) to 0.35 eV (with a step of 0.05 eV). The ordering 
temperatures as a function of the band overlaps are shown 
in the insets. The lower T c for J=0.1 eV and A=0.05eV 
correcpond to the second phase transition. 



lap of about 0.3 eV. The selfconsistent GW calculation 12 
finds CaB 6 to be a semiconductor with a band gap of the 
order of 1 eV, as does the more conventional pseudopo- 
tential GW calculation*^ The GW method is known to 
provide good bandgaps for many semiconductors where 
the LDA gaps are strongly underestimated or even van- 
ishing. The effect of GW self-energy corrections on the 
EuBg bandstructure is expected to be similar, thus re- 
ducing the band overlap and possibly opening a band 
gap by shifting the conduction band upwards. 

An important difference between EuBg and CaBg is the 
presence of 4/ orbitals. We have shown that the exchange 
interaction with the band electrons is opposite in sign 
for the valence and conduction bands, so magnetic order 
can significantly increase the overlap of the majority spin 
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band while it has opposite effect for the minority spin. 
Although the effective exchange with the valence band 
depends on the position of the 4/ states and thus varies 
somewhat with the parameter U, general experience with 
rare earth systems as well as optical measurements on 
EuBgS indicate that U of 7 eV (placing the Eu 4/ 1 eV 
below the Fermi level) is realistic. 

Based on these arguments it is plausible that the realis- 
tic ground state picture of stochiometric EuBg is that of a 
half metallic semimetal. This unprecedented band struc- 
ture would result from a ferromagnetic GW calculation 
in which the band shifts (relative to the static LDA+U) 
are large open a gap in the minority spin channel, but 
not so large as to open a gap in the majority spin chan- 
nel. Such a scenario has some experimental support, dc 
Haas-van Alphen data provide clear indication of Fermi 
surfaces, but only two pockets are seen^iSS This num- 
ber of sheets is contrary to the expected four pockets 
suggested by LDA (and LDA+U) band structure, but 
two sheets are exactly what is expected of a half metallic 
semimetal. 

The observed behavior of the resistivity is also consis- 
tent with such a picture. In the magnetically ordered 
state, metallic conduction takes place in the majority 
spin channel. Just above the magnetic ordering tem- 
perature the system can be viewed as consisting of dis- 
ordered magnetic domains (due to short range ferromag- 
netic correlations), and the increase of resistivity upon 
disordering is due to mismatch of the conducting spin 
channels between these domains (like the "intergrain tun- 
neling mechanism" of giant magnetoresistance materials; 
see e.g. Ref. 12911 . Increasing the temperature further 
leads to breakdown of the short range order, the para- 
magnetic bandstructure becomes increasingly appropri- 
ate, and the inter-domain magnetoresistance effect dis- 
appears, as observed. 

Angle-resolved photoemission datai^ do not seem to fit 
well into this picture, but further photoemission studies 
in the ordered phase, and identifying the position and 
influence of the / states, seems to be necessary to clarify 
several remaining questions. The very large number of 
electron carriers that the photoemission data, if assumed 
to be representative of the bulk, does not fit so well with 
data that suggest rather clean single crystals. If the car- 
riers are due to unbalanced surface charge, then the ob- 
served bands are not representative of the bulk. Finally, 
the broken inter-octahedron B-B bonds should give rise 
to surface states (or bands), and the implication of the 
photoemission data will never be unambiguous until the 
surface electronic structure is identified and understood. 

Finally we discuss briefly our mean field treatment of 
the two band Kondo lattice model hamiltonian. We ex- 
pect that the least reliable quantity obtained from the 
electronic structure calculation is the band overlap, and 
we have treated it as a parameter in the mean field study. 
The two band hamiltonian even in mean field approxima- 
tion leads to a strongly temperature dependent coupling 



between the local moments, which is reflected by unusual 
magnetization dependences distinct from the canonical 
behavior of the Heisenberg hamiltonian in mean field. 
The solutions for small positive band overlap exhibit a 
slow approach to saturation as observed^ as well as an 
increase of the ordering temperature with pressure^ i.e. 
with increasing band overlap. The ordering temperatures 
obtained with the ab inito value of the exchange param- 
eter J are, however, too low compared to experiment. 
We mention several possible reasons for this discrepancy. 
First, we have completely neglected the k-dependence of 
the exchange parameter J. Second, equation (8) shows 
that the value of J v is quite sensitive to the position 
of the occupied 4/ levels, which we know only approx- 
imately. And finally, some mechanisms of effective ex- 
change between the local moments arising from the An- 
derson model are lost when transformed to the Kondo 
model (e.g. superexchange) . These mechanisms might 
be of importance in the A < case leading to removal of 
the first order transition. The mean field study does serve 
to demonstrate the potential usefulness of the two-band 
hamiltonian and to stimulate further studies by more ad- 
vanced techniques, such as the Green's function approach 
of Nolting et alM 



V. CONCLUSIONS 

Using the LDA+U approach we have shown that treat- 
ing the Eu 4/ states within the same framework as 
the rest of the itinerant electrons has important conse- 
quences. In particular it leads to Kondo coupling be- 
tween local moments and valence electrons, but anti- 
Kondo coupling to conduction electrons. We have iden- 
tified and described in detail the origin of this coupling 
as well as the origin of the k-dependence of the anti- 
Kondo (ferromagnetic) coupling to the conduction band. 
Based on our electronic structure analysis we suggest the 
description of EuBg in terms of two band Kondo/anti- 
Kondo lattice model, and shown that a half metallic 
semimetal results at the mean field level of description. 
We have obtained the parameters of the corresponding 
hamiltonian, which will allow more extensive material- 
specific treatments in the future, and demonstrated the 
effects on magnetic ordering arising from an exchange 
controlled band overlap. The picture we present seems 
consistent with observed Fermi surfaces and transport 
properties. 
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